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^ ■ 1 Introduction 

o 

D ' 1.1 Motivation and definitions 

Q" 

The topic of spacing distributions in random matrix ensembles is almost as old as the introduc- 
ed ' tion of random matrix theory into nuclear physics. Both events can be traced back to Wigner 
in the mid 1950's |37l I38j . Thus Wigner introduced the model of a large real symmetric random 
matrix, in which the upper triangular elements are independently distributed with zero mean 
\Q \ and constant variance, for purposes of reproducing the statistical properties of the highly ex- 
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cited energy levels of heavy nuclei. This was motivated by the gathering of experimental data 
on the spectrum of isotopes such as 238 U at energy levels beyond neutron threshold. Wigner 
hypothesized that the statistical properties of the highly excited states of complex nuclei would 
be the same as those of the eigenvalues of large random real symmetric matrices. For the 
random matrix model to be of use at a quantitative level, it was necessary to deduce analytic 
forms of statistics of the eigenvalues which could be compared against statistics determined from 
experimental data. 

What are natural statistics for a sequence of energy levels, and can these statistics be com- 
puted for the random matrix model? Regarding the first question, let us think of the sequence as 
a point process on the line, and suppose for simplicity that the density of points is uniform and 
has been normalized to unity. For any point process in one dimension a fundamental quantity 
is the probability density function for the event that given there is a point at the origin, there 
is a point in the interval [s, s + ds], and further there are n points somewhere in between these 
points and thus in the interval (0, s). Let us denote the probability density function by p(n; s). 
In the language of energy levels, this is the spacing distribution between levels n apart. 

Another fundamental statistical quantity is the fc-point distribution function pr k \{x\, . . . , x^j. 
This can be defined recursively, starting with p^{x), by the requirement that 

P(k){ x U-- ■ : x k)l P(k—l) O^li • • ■ , x k-l) (1-1) 

is equal to the density of points at x^, given there are points at x\, . . . , x^-i- One sees from the 
definitions that 

P(2)(M _f> , x (12) 
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From empirical data of a long energy level sequence, the quantity p(n; s) for small values of n 
at least is readily estimated (the statistical uncertainty gets worse as n increases). Use of (jl.2|) 
then allows for an estimation of /?(2)(0; s). 

We thus seek the theoretical determination of p(n; s) for matrix ensembles. 

1.2 Spacing between primes 

Before taking up the problem of determining p{n\ s) for matrix ensembles, which is the theme of 
these lectures, let us digress a little and follow the line of introduction to spacing distributions 
given by Porter in the review he wrote as part of the book jSJ, which collected together the 
major papers written in the field up to 1965. Porter's introduction is particularly relevant to the 
theme of the present school because it uses the prime numbers as an example of a deterministic 
sequence which, like energy levels of heavy nuclei, nevertheless exhibit pronounced stochastic 
features. 

It turns out the spacing distributions between primes relate to perhaps the simplest example 
of a point process. This is when the probability that there is a point in the interval [s, s + ds] 
is equal to ds, independent of the location of the other points. This generates the so called 
Poisson process with unit density, or in the language of statistical mechanics, a perfect gas. By 
definition of the process the ratio 1)1.1(1 is unity for all k and thus 

P(k)(xi,...,x k ) = 1. (1.3) 

To compute p(n; s), we think of the Poisson process as the N — * oo limit of a process in which 
each unit interval on the line is broken up into N equal sub-intervals, with the probability of 
there being a particle in any one of the subintervals equal to 1 /N. Thus 

1 /q/V\ t n 

p(s; n) = lim (1 - —) sN - n N- n ( ) = ^e" 8 . (1.4) 
n^oo N \ n J n! 

In the first equality of (|1.4|) . the first factor is the probability that sN — n subintervals do 
not contain a particle, the second factor is the probability that n subintervals do contain a 
particle, while the final factor is the number of ways of choosing n occupied sites amongst sN 
sites in total. The probability density in the final equality of (|l,4j) is the Poisson distribution. 
Substituting (|1.4|) in (|1.2|) gives p(2)(0,x) = 1, as required by p.3jl . 

The distribution (|1.4j) ties in with prime numbers through Kramer's model (see the lectures by 
Heath-Brown in the present volume). In this approximation, statistically the primes are regarded 
as forming a Poisson process on the positive integer lattice. The probability of occupation of 
the iVth site is taken to equal 1 / log N, so as to be consistent with the prime number theorem. 
Kramer's model predicts that as an approximation 

p (A,) (n;s) = i" s , s = i/log7V (1.5) 
n! 

where p( N \n; s) refers to the probability that for primes p in the neighbourhood of a prime N, 
there is a prime at p + i, and furthermore there are exactly n primes between p and p + t. 
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Figure 1: Distribution of the spacing t between primes (leftmost graph) and the spacing t 
between every second prime for 2, 000 consecutive primes starting with N = 10 9 + 7. The 
distributions are given in units of s = tj log N. The smooth curves are the Poisson distributions 
p(0; s) = e~ s and p(l; s) = se~ s . 

To compare the prediction H1.5|) against empirical data, we choose a value of N, say 10 9 , 
and for the subsequent M primes (say M = 2, 000) record the distance to the following prime 
(in relation to p( N \l;s)) and the distance to the second biggest prime after that (in relation 
to p( N '(s; 1)). We form a histogram, with the scale on the horizontal axis measured in units of 
s = t/logN, where t is the actual spacing. The natural units for t are multiples of 2, and this 
provides a width for the bars of the histogram. We see from Figure ^ that the general trend of 
the histograms do indeed follow the respective Poisson distributions. 

1.3 Empirical determination of spacing distributions for matrix ensembles 

Wigner's interest was in the statistical properties of the eigenvalues of large real symmetric 
random matrices. More particularly, he sought the statistical properties of the eigenvalues in 
what may be termed the bulk of the spectrum (as opposed to the edge of the spectrum [S]). The 
eigenvalues in this region are characterized by having a uniform density, which after rescaling 
(referred to as 'unfolding') may be taken as unity. In distinction to the situation with the 
sequence of primes, for random matrices it is not necessary to study the statistical properties of 
a large sequence of (unfolded) eigenvalues from a single matrix. Rather the spacing distributions 
with respect to the middle eigenvalue (this is the eigenvalue most in the bulk) in multiple samples 
from the class of random matrices in question can be listed, and then this list used to create 
a histogram. Moreover, to approximate large matrix size behaviour, it is only necessary to 
consider quite small matrix sizes, say 13 x 13. 

In Figure[2]we have plotted the empirical determination of p(0; s) and p(l; s) obtained from 
lists of eigenvalue spacings for realizations of the so called GOE (Gaussian orthogonal ensemble) 
eigenvalue distribution. As we know from the lectures of Fyodorov in this volume, the GOE 
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Figure 2: Plot of the distribution of the unfolded spacing between the 6th and 7th, and 7th and 
8th eigenvalues (pooled together) for 2,000 samples from the 13 x 13 GUE eigenvalue distribution. 
The smooth curve is the Wigner surmise (|2.12jl . The rightmost graph is the distribution between 
the 6th and 8th eigenvalues in the same setting, while in the smooth curve in this case is 
(l/2)p 4 (0; s/2) with p 4 given by (f2~T3|) . 



consists of real symmetric random matrices, with each diagonal element chosen from the normal 
distribution N[0, 1], and each (strictly) upper triangular element chosen from the normal distri- 
bution N[0, l/v2]- For such matrices, it is well known that to leading order in the matrix rank 
N, the eigenvalue density is given by the Wigner semi-circle law 



P(i)O) = 



2N / x 2 



Multiplying the eigenvalues at point x by this factor allows us to unfold the sequence giving a 
mean eigenvalue spacing of unity 

A less well known, and much more recent result relating to GOE matrices is that their 
spectrum can be realized without having to diagonalize a matrix jlj (see also |13j). Thus one 
has that the roots of the random polynomial P/v(A), defined recursively by the stochastic three 
term recurrence 

P k (X) = (A - a k )P k ^{\) - &L 1 P fc - 2 (A) (1.6) 

where 

a k ~ N[0, 1], b\ ~ Gamma[Jfe/2,l], 

have the same distribution as the eigenvalues of GOE matrices (the notation Gammafs, a] de- 
notes the gamma distribution with density proportional to x s ~ 1 e~" x / a ). Generating such polyno- 
mials and finding their zeros then provides us with a sequence distributed as for GOE eigenvalues, 
from which we have determined p(0; s) and p(l; s). 
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2 Eigenvalue product formulas for gap probabilities 

2.1 Theory relating to p(n; s) 

Consider a point process consisting of a total of N points. Let the joint probability density 
function of the N points be denoted p(xi, ... ,xn). A quantity closely related to the spacing 
distribution p(0; s) is the gap probability 

E h (0;s):= lim / dx\--- I dx^ p{&N x u ■■■ > &n%n) (2-1) 

TV— oo Jj Jj 

where I = (—00,00) — (— s/2, s/2) and a^v is the leading large N form of the local density at 
the origin (and thus the unfolding factor). Thus it is easy to see that 

p(0;s) = -^E h - lk (0;s). (2.2) 

More generally we can define 

E hulk (n; s) := lim ( )a N / dx%--- / dx n I dx n+ i--- / dxjsrpia^xi, . . . , oat^at)- 
N^oo \ n J J_ s/2 J_ g/2 J j J j 

(2.3) 

These quantities can be calculated from the generating function 



/oo poo N 

dx x ■ ■ ■ / dx N - CX I s/2iS/2 ))p(«iV^i> • • • , <*nxn), (2.4) 
-00 J — 00 

where Xj = 1 f° r x ® 6 J and Xj = otherwise, according to the formula 

^ lk (n; S ) = i- | L— ^ lk ( S ;e)| 5=i . (2.5) 

It follows from the definitions that 

d 2 

p(n- s) = —E hnlk (n; s) + 2p(n - 1; s) - p(n - 2; s), (2.6) 
ds z 

or equivalently 

pfa s ) = j~2 D ra - j + !)^ buIk (i; *)■ ( 2 - 7 ) 

S i=o 

Hence knowledge of {E h (j; s)}j=o,...,n is sufficient for the calculation of p(n;s). 

It is possible to relate Q2.4J1 to the fc-point distribution functions. In the finite system the 
latter are given by 

P(fc) (^1) • • • > x k) = , N J dx k+1 ---J dx N p(x 1} ... ,x N ). (2.8) 

With 

bulk/ \ t b (N) 1 \ 

pfy (xi,...,x k ) := hm a N p^ (a N xi, . . . ,a N x k ), 
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by expanding 1)2.4(1 in a power series in £ and making use of (|2.8|) we see that 

\k rs/2 i-a/2 



k— 1 ' J— s/2 J— s/2 



(2. 



For the limiting process to be rigorously justified, because [—s/2, s/2] is a compact interval, it 
is sufficient that pj?^ lk (xi, . . . , x k ) be bounded by M k for some M > 0. 

With these basic formulas established, we will now proceed to survey some of the main 
results relating to spacing distributions in the bulk of the various matrix ensembles (orthogonal, 
unitary and symplectic symmetry classes). 

2.2 Wigner surmise 

For the Poisson process we have seen that p(0; s) = e~ s . Thus in this case the spacing distribution 
is actually maximum at zero separation between the points. The opposite feature is expected 
for p(0; s) in relation to the eigenvalues of random real symmetric matrices, as can be seen by 
examining the 2x2 case of matrices of the form 



A 



a b 
b c 



This matrix is diagonalized by the decomposition A = i?diag[A+, \-]R T where 



R 



COS I 

sin ( 



— sin ( 

cos 6 



Expressing a, b, c in terms of A+, A_, 8 it is simple to show 

dadbdc = |A+ — \-\d\+d\-d8. 



(2.10) 



Thus for small separation s := |A+ — A_| the probability density function vanishes linearly. 

Let fj,(s) denote the small s behaviour of p(0;s). We have seen that for the Poisson pro- 
cess /j,(s) = 1, while for the bulk eigenvalues of real symmetric matrices /x(s) oc s. Wigner 
hypothesized that as with the Poisson process, p(0; s) for the bulk eigenvalues of random 
real symmetric matrices could be deduced from the ansatz 

p(0; s) = cin(s) exp ( - c 2 J n(t) dtj (2.11) 

where the constants c\ and c 2 are determined by the normalization requirements 

/■OO /'OO 

/ p(0; s) ds = 1, / sp(0;s)ds = l 
Jo Jo 

(the second of these says that the mean spacing is unity). Thus one arrives at the so called 
Wigner surmise 



p(0;s) 



*-se-^l* 



(2.12) 
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for the spacing distribution of the bulk eigenvalues of random real symmetric matrices. 

The ansatz (|2.11j) does not apply if instead of real symmetric matrices one considers complex 
Hermitian matrices, or Hermitian matrices with real quaternion elements. Examining the 2x2 
case (see the introductory article by Porter in |31j ) one sees that in the analogue of (|2.1UI) . the 
factor |A + — A_| should be replaced by |A + — A_|^ with (3 = 2 (complex elements) or (3 = 4 
(real quaternion elements). Choosing the elements to be appropriate Gaussians, one can reclaim 
(|2.12() and furthermore obtain 

p 2 (0; s) = ^e-^l\ p 4 (0; s) = ^e^ 2 ^. (2.13) 

7T Z 3 D 7T J 

as approximations to the spacing distributions in the cases (3 = 2 and (3 = 4 respectively. 



2.3 Predholm determinant evaluations 



A unitary invariant matrix ensemble of N x N random complex Hermitian matrices has as its 
eigenvalue probability density function 

N 



1=1 Kj<k<N 



(2.14) 



which we will denote by UEtv(#). We know (see the lectures by Fyodorov in this volume) that 
the fe-point distribution function can be expressed in terms of the monic orthogonal polynomials 
{Pk(x)}k=o,i,... associated with the weight function w 2 (x), 



w 2 (x)pj(x)p k (x)dx = hj5j tk - 



Thus with 



7V-1 



K N ( X ,y) = (w 2 {x)w 2 (y)) 1 / 2 £ PJ&M 



k=0 



(w 2 (x)w 2 {yj) 



i/ 2 Pn{x)pn-i{v) -pN(y)PN-i(x) 
x-y 



we have 



p^\xi, . . . ,x k ) = det K N (xj,xi) 



j,i=i,...,k 

This structure is significant for the evaluation of the generating function 



N 



E N , 2 (J;Z;w 2 ) :=m(l-£ x «) 



l=i 



XJE N (g) 



(2.15) 
(2.16) 

(2.17) 



(the subscript 2 on En,2 indicates the exponent in (|2.14|) ). Expanding (|2.17|) in a power series 
analogous to ()2.9|) we obtain 



A' 



k=l 



dx-\ 



dxu det 



K N (xj,xi) 



j,l=i,...,k 



(2.18) 
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where use has been made of (|2.16|) . The sum in (|2,18|) occurs in the theory of Fredholm integral 
equations [HEl; and is in fact an expansion of the determinant of an integral operator, 



E N>2 {J-i-w 2 ) = det(l-eFO) 
where Kj is the integral operator on the interval J with kernel K^(x,y), 



(2.19) 



K N [f](x) 



KN{x,y)f{y)dy. 



It is well known that in the bulk scaling limit, independent of the precise functional form of 

w 2 (x), 

din 7r( nr — 7/1 , 

\x,y) (2.20) 



. , sin7r(x — y) . 
lim a N K N {a N x,a N y) = ^ ^ =: K bull > 

N^oo TT[X — y) 



for a suitable scale factor a at- Thus we have 

E% nlk (J;£) = det(l - e^j ulk ) 



(2.21) 



where Kj n]k is the integral operator on the interval J with kernel (|2.20|) (the so called sine kernel). 
This is a practical formula for the computation of Evp^ if we can compute the eigenvalues 
{/^■}j=o,i,... of Kj ulk , since we have 



^2 bulk («/;0 = II( 1 -^ 

j=0 



(2.22) 



In fact for J = (— s, s) the eigenvalues can be computed ^H] by relating K^ lk s ^ to a differential 
operator which has the prolate spheroidal functions as its eigenfunctions, and using previously 
computed properties of this eigensystem. 

Wigner's interest was not in complex Hermitian random matrices, but rather real symmetric 
random matrices. Orthogonally invariant ensembles of the latter have an eigenvalue probability 
density function of the form 

1 N 

^Il^ifa) II \ x k-Xj\, (2.23) 



1=1 



l<j<k<N 



to be denoted OEjv(wi)- For such matrix ensembles, the Appoint distribution function can be 
written as a quaternion determinant (or equivalently Pfafnan) with an underlying 2x2 matrix 
kernel (see e.g. jSJ Ch. 5]). From this it is possible to show that 



^ibull; 



(J;0 



det(l - iK\ 



bullo 



(2.24) 



where K b u J k is the integral operator on J with matrix kernel 

smir(x-y) 1 r (l_!,) sinf 



R bu\k 



tt(x — y) 7T 
d sin7r(x — y) 
dx tt(x — y) 



/ dt - -sgn(x - y) 

sin7r(x — y) 
tt(x - y) 



(2.25) 
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However, unlike the result (|2.21|) . this form has not been put to any practical use. 

Instead, as discovered by Mehta [2E], a tractable formula results from the scaling limit of an 
inter-relationship between the generating function of an orthogonal symmetry gap probability 
and a unitary symmetry gap probability. The inter-relationship states 

E 2N ,i{(-t, t)\ £; e" xV2 )| 5=1 = E Nt2 ((0, t 2 )- £; y -i/2 e -» Xl , >0 )| (2.26) 

and in the scaling limit leads to the result 

E^ lk ((s,s);0 ^ = det(l - K^+) (2.27) 



where Kr^j\ is the integral operator on (—s, s) with kernel 



1 /sin-7r(x — y) sin7r(x + y)\ 

+ — 7^ — rN> ( 2 - 28 ) 



2 V ir(x — y) n(x + y) 

which we recognize as the even part of the sine kernel (12.20(1 . (For future reference we define 
~s analogously, except that the kernel consists of the difference of the two terms in (|2.28|) , 
or equivalently the odd part of the sine kernel (|2.20|) .) Because the eigenvalues /j,2j of the integral 
operator on (— s, s) with kernel (|2.2U|) correspond to even eigenfunctions, while the eigenvalues 
fj>2j+i correspond to odd eigenfunctions, we have that 

oo 

Sf nIk ((-« ) -);OL =1 = II(l-Mai)- (2-29) 

5 1=0 

Gaudin ^H] used this formula, together with (|2.2|) . to tabulate p\ nlk (0; s) and so test the accuracy 
of the Wigner surmise (J2.12|) . In fact this confirmed the remarkable precision of the latter, with 
the discrepancy between it and the exact value no worse than a few percent. 

The case of Hermitian matrices with real quaternion elements and having a symplectic sym- 
metry remains. The eigenvalue p.d.f. of the independent eigenvalues (the spectrum is doubly 
degenerate) is then 

1 N 

^n^(^) n (^-^) 4 > (2-30) 

1=1 l<j<k<N 

which we denote by SEn(w4 : ). The computation of the corresponding bulk gap probability relies 
on further inter-relationships between matrix ensembles with different underlying symmetries. 
These apply to the eigenvalue probability density function for Dyson's circular ensembles, 

l<j<k<N 

where = 1,2 or 4 according to the underlying symmetry being orthogonal, unitary or sym- 
plectic respectively. The corresponding matrix ensembles are referred to as the COEjv, CUEjv 
and CSEjy in order. In the N — > oo scaling limit these ensembles correspond with the bulk of 
the ensembles OEjv(itfi), \JEn(w2) and SEa^u^) respectively. 
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The first of the required inter-relationships was formulated by Dyson [3] and proved by 
Gunson JH]- It states that 

alt(COE iV U COEtv) = CUE^v (2.31) 

where the operation COE^r U COEtv refers to the superposition of two independent realizations 
of the COEtv and alt refers to the operation of observing only every second member of the 
sequence. The second of the required inter-relationships is due to Dyson and Mehta [7]. It 
states that 

altCOE 2 iv = CSE^v. (2.32) 

(For generalizations of (|2.31l) and (|2.32|) to the ensembles OEn(w\), \JEn(w2) and SEtv(^4) 
with particular w\, W2 and see [12].) Using Q2.31fl and (|2.32|) together one can deduce that 
in the scaled limit 

i?r k (0; (s/2, s/2)) = \ (EF* (0; (s, a)) + ^^^) , (2-33) 
which upon using (|2.22f) and (|2.29f) reads 

^ 00 00 
< lk (0; (s/2, s/2)) = - ( ]J(1 " hi) + H(l - Aa+0) . (2.34) 

1=0 1=0 

Another consequence of (|2.32j) is that 

p 4 (0;a) = 2 Pl (l;2a). (2.35) 

It is this relationship, used together with the approximation for pi(0; s) in (|2.13j) . which is used 
to approximate p(l; s) smooth curve in Figure [2 

In summary, as a consequence of the pioneering work of Mehta, Gaudin and Dyson, com- 
putable formula in terms of the eigenvalues of the integral operator on (— s,s) with the sine 
kernel (|2.20|) were obtained for 

E% nl H(s, s);0, £i ulk (0; (s, s)), £ 4 bulk (0; (s/2, s/2)). 

3 Painleve transcendent evaluations 
3.1 The results of Jimbo et al. 

An explicit connection between the multiple interval gap probability 

and integrable systems theory — specifically the theory of isomondromic deformations of linear 
differential equations — was made by Jimbo, Miwa, Mori and Sato in 1980. Here the endpoints 
ai, . . . , 02 P of the gap free intervals become dynamical time like variables, inducing flows which 
turn out to be integrable. 
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As part of this study the quantity 

oo 

E™ k ((-s,s);0 = det(l - ZK*% a) ) = - frj) (3.1) 

3=0 

was expressed in terms of the solution of a nonlinear equation. In fact knowledge of ()3.1|) is 
sufficient to calculate the products appearing in (|2.29f) and (|2.34|) . Thus with 

oo oo 

D + (s;0 := H(l - & 2j ), D_( 8 ;0 := " Z&j+l) 

3=0 j=0 

Gaudin (see has shown 



1 1 f s / d 2 

log J D±( S ;6 = 2 lo g^2 bulk ((-^^);6±2 J Q \J-^2^gE^H(-x,x);0dx. (3.2) 

The result of [23] is that 

E^ k ((-s,s);0 =exp H du (3.3) 

J o u 

where cr(u; £) satisfies the nonlinear differential equation 

{uu"f + 4(ua' - a){ua' -a + (a') 2 ) = (3.4) 
subject to the boundary condition 

U-+0+ 7T 

In fact the equation 1)3. 4 J) is an example of the so called a form of a Painleve V equation. In 
view of this it is appropriate to give some background into the Painleve theory, following |21j . 
First we remard that the Painleve differential equations are second order nonlinear equations 
isolated as part of the study of Painleve and his students into the moveable singularities of 
the solution of such equations. Earlier Fuchs and Poincare had studied first order differential 
equations of the form 

P(y',y,t) = (3.5) 

where P is a polynomial in y 1 , y with coefficients meromorphic in t. In contrast to linear differ- 
ential equations, nonlinear equations have the property that the position of the singularities of 
the solution will depend in general on the initial condition. The singularities are then said to 
be moveable. For example 

y 2 (3.6) 



dt 

has the general solution y = l/(c — t), where c determines the initial condition, and so exhibits 
a moveable first order pole. The nonlinear equation 

dy = l 
V dt 2 
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has the general solution y = (t — c) 1 / 2 , which exhibits a moveable branch point (essential 
singularity) . Fuchs and Poincare sought to classify all equations of the form (|3.5|) which are free 
of moveable essential singularities. They were able to show that up to an analytic change of 
variables, or fractional linear transformation, the only such equations with this property were 
the differential equation of the Weierstrass "P-function, 



% 3 - <? 2 y -S3, (3.7) 



dtJ 

or the Riccati equation 



^ = a(t)y 2 + b(t)y + c(t) (3.8) 

where a, b, c are analytic in t (note that (|3.6j) is of the latter form). 

Painleve then took up the same problem as that addressed by Fuchs and Poincare, but now 
with respect to second order differential equations of the form 

y" = R(y',y,t) 

where R is a rational function in all arguments. It was found that the only equations of this form 
and with no moveable essential singularities were either reducible to Q3.7|) or l|3.8|) . reducible to 
a linear differential equation, or were one of six new nonlinear differential equations, now known 
as the Painleve equations. As an explicit example of the latter, we note the Painleve V equation 
reads 

/ = (1 + > ) [v <f - -,J + fclll! U + Z) + 22 + 5 Jik±21 (3.9) 
\2y 1 — yJ x ar V yJ x y — 1 

where a,/3,7 are parameters. 

An immediate question is to how (|3.9j) relates to (|3.4| . For this one must develop a Hamil- 
tonian theory of the Painleve equations. The idea is to present a Hamiltonian H = H(p, q, t; v), 
where the components of v are parameters, such that after eliminating p in the Hamilton equa- 
tions 

. dH . dH . 

" = w W (3 ' 10) 

q' and p' denoting derivatives with respect to t, the equation in q is the appropriate Painleve 
equation (in p.lUj) the role of p and q is interchanged relative to their usual meaning of position 
and momentum in physics; here we are following the convention of Okamoto. Malmquist 
was the first to present such Hamiltonians, although his motivation was not to further the 
development of the Painleve theory itself. This was left to Okamoto in a later era, and it is 
aspects of his theory we will briefly present here. 

The Hamiltonian for the PV equation as presented by Okamoto 29 is 

tH v = q(q-l) 2 p 2 - {(V! - v 2 ){q - I) 2 - 2{v x + v 2 )q{q - 1) + tq}p 

+ (v3-v 2 )(v4-v 2 )(q-l), (3.11) 

where the parameters are constrained by v± + v 2 + v 3 + V4 = and are further related to those 
in (|3.9I) according to 

a = ^(v 3 - v 4 ) 2 , (3 = ~^(vi - v 2 ) 2 , 7 = t?i + 2v 2 - 1, 5 = 
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It turns out that, as a consequence of the Hamilton equations (j3.1U|) . tHy itself satisfies a 
nonlinear differential equation. It is this differential equation which relates to (j3.4|h Okamoto 
made use of this equation for the symmetry it exhibits in the parameters v\, . . . , U4. 

The equation in question, which is fairly straightforward to derive, is presented for the so 
called auxilary Hamiltonian 

h v {t) = tHy + (v 3 - V 2 )(v 4 - v 2 ) - v 2 t - 2v\. 

Okamoto showed 

4 

{th'yf - (h V - th'y + 2{h'y) 2 ) 2 + 4 J] (h'y + V k ) = 0. 

fc=l 

Setting 

a v {t) = hy(t) + v 2 t + 2v 2 , Vj-i = Vj -v 2 (j = 1, . . . , 4) 
in this one obtains the so called Jimbo-Miwa-Okamoto cr-form of the Painleve V equation 

{ta'yf - (ay - ta' v + 2{a' v ) 2 + (u + Vl + u 2 + ^Wy) 2 

+4(u + a'yXu! + a' v )(u 2 + a' v )(u 3 + a' v ) = (3.12) 

(Jimbo and Miwa arrived at (|3.12j) in their study of isomonodromic deformations of linear 
differential equations). We note that (|3.4|) is an example of this equation with 

vo = V\ = v 2 = f3 = 0, 1 1 ► —2iu. 
3.2 Unveiling more structure 

The result of Jimbo et al. relates to the Fredholm determinant of the integral operator with the 
sine kernel. What is special about the sine kernel that relates it to integrable systems theory? 
This question was answered by Its, Izergin, Korepin and Slanov |2U| who exhibited integrability 
features of all kernels of the Christoffel-Darboux type (recall (|2. 15f) in relation to the latter 
terminology) 

tK(x, y) = ^HM ; (3.13) 
x-y 

the sine kernel begin the special case 

4>{x) = v^sinx, ip(y) = y^cosy. (3-14) 

One of their key results related to the form of the kernel R(x, y) for the so called resolvent 
operator 

Rj :=iKj(l-iKj)-\ 

With 

Q(x) := (1 - ZKj)- 1 ^ P(x):={l-iKj)- 1 ^ (3.15) 
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they showed 

^^.^W-i-MW, (3 . 16) 

x - y 

The significance of the resolvent kernel is evident from the general formula 

Ai ogdet (l-£if (ai , 02) ) = (-iy- l R( aj , aj ) (j = 1,2). (3.17) 
To derive this formula, one notes that 

logdet(l-^ (ai)02) ) = Tr log(l-^ (oiiaa) ) 

(s) 



Thus 



tog(l-^if(a;,x)xj a ; )tta) )dx 



^-logdet(l-£K (ai , a2) ) = (-iy- 1 (l-^(a i ,a i ))- 1 ^(a i ,, 



as required. 

According to (|3.16j) 



R{aj,aj) = -Q{x)P'{x) + P(x)Q'(x) , (3.18) 



so we see from (j3.17|) that the Fredholm determinant is determined by the quantities (|3.15|) and 
their derivatives evaluated at the endpoints of the interval. Indeed a close examination of the 
workings of [23] , undertaken by Mehta [27] , Dyson [Hj and Tracy and Widom [32] , revealed that 
the former study indeed proceeds via the equations (|3.17|) and 1)3.18(1 . and in fact a(t) in ()3.3(1 
is related to the resolvent kernel evaluated at an endpoint by a(t) = —tR(t/2,t/2). Moreover 
it was realized that like ((3.16J) there are other equations contained in the working of |23j which 
apply to all kernels of the form 1)3.13(1 . However it was also clear that other equations used in 
|23j were specific to the form of <f> and ip in 1)3 .14(1 . 

Tracy and Widom were able to identify these latter properties, which are that 4> an d ip are 
related by the coupled first order differential equations 

m{x)4>'{x) = A{x)4>{x) + B(x)i(j(x) 

m{x)ip'(x) = -C(x)<p(x) - A(x)ip(x) (3.19) 

where m, A, B, C are polynomials. This structure allows the so called universal equations (in- 
dependent of the specific form of (|3.13|0 such as (|3.18|) to be supplemented by a number of case 
specific equations. For some choices of (ft and i/j in addition to that corresponding to sine kernel, 
the resulting system of equations closes. Examples relevant to spacing distributions at the soft 
and hard edge of matrix ensembles with unitary symmetry are 

<j){x) = \/lki(x), ip(x) = 4>'{x), <j)(x) = \flJ a {yJ~x), ip(x) = x(j)'(x). 

In both these cases it was possible to obtain an evaluation of the generating function for the 
corresponding gap probability in a form analogous to (|3.3|) 33^ . 



14 



We will make note of the hard edge result because it, by virtue of Mehta's inter-relationship 
(|2.2fi|) , relates to the gap probability in the bulk in the case of an underlying orthogonal symme- 
try. First, we define the hard edge gap probability in the case of an underlying unitary symmetry 
as the scaled limit of the ensemble (j2.14|) with W2(x) = x a e~ x Xx>o- Explicitly 

£ 2 hard ((0,s);£) = lim E 2 ((0,^-y,^x a e- x Xx >o). (3.20) 



£ 2 hard ((0, S );£) = det(l-£K» (3-21) 



It was shown in [2j that 

L (CM) 

where K^jl is the integral operator on (0, s) with kernel 

K ^ d( v _ UxW)yV*J' a {yW) - x^J' a (x^)J a (y^) 
[X ' y) ~ 2{x-y) 
As part of the study [33] the Fredholm determinant (|3.21|) was given the evaluation 

f s dt 

E^ d ((0,s);0=exp u(i;a;£)- (3.22) 
Jo 1 

where u satisfies the differential equation 

itu") 2 - a 2 (u') 2 - u'(4u' + l)(u - tv!) = (3.23) 

subject to the boundary condition 

u(t; a; ~ -£tK hard (t, t). 
t^o+ 

The equation (|3.23|) is a special case of the a form of the Painleve III' system |30| . 
It follows from (t2~2ll. (l3~2771) and (l3~22l that |XD| 

/*O s ) 2 A+ 

E™"(0; (-8,8)) = £ 2 hard (0; (0,vrV)) = exp | u(t;a;0 j __ i/2 . (3.24) 

This is an alternative Painleve transcendent evaluation to that implied by (|2.28|) , (|3,2j) and l)3.3j) . 
Similarly, by noting that 

l/sin(x — y) sin(x + y) ' 
a=i/2 2 V x — y x + y 



2^K h ™ d {x\y 2 ) 



we see from lf2~52jl . (l3~2TT) and (l3~22l that H 

E^ lk (0;(-s/2,s/2)) 

1 / /'( 7rs ) 2 dt ff 18 ) 2 dt \ 

= 2 v exp j a; ^ T p + exp 7 a; ^ T *=_\ /2 ) ' (3 ' 25) 

In summary, the Fredholm determinants in the expressions for the bulk gap probabilities 
can each be written in terms of Painleve transcendents. From a practical viewpoint these 
expressions are particularly well suited for generating power series expansions, and also allow for 
a numerical tabulation of each of £J? ulk (0; (-a, a)), E\ ulk (0; (-a, s)) and E% ulk (0; (-a, s)), as well 
as i?2 Ulk (^; (— s,s)) for n > 1. For the latter quantity, according to (|2.5[) we must differentiate 
E\ (( — s ; s )jC) with respect to £ then set £ = 1. Doing this in (|3.4|) gives a coupled system 
of differential equations for d 3 a{u; £)/5£ J |^ = i (j = 0, . . . ,n) which is only numerically stable for 
small values of n. 
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3.3 Distribution of bulk right or left nearest neighbour spacings 

The spacing distribution refers to the distribution of the distance between consecutive points 
as we move along the line left to right. Another simple to measure statistic of this type is the 
distribution of the smallest of the left neighbour spacing and right neighbour spacing for each 
point. Let us denote this by Pp' n '(s) (the superscript n.n. stands for nearest neighbour, while 
the subscript indicates the symmetry class). Let £^' n '(0; (— s,s)) denote the probability that 
about a fixed eigenvalue at the origin, there is no eigenvalue at distance s either side. Analogous 
to (|2.2j) it is easy to see that 

^' n ' (s) = "^' n ' (0;( " S ' s)) - (3 - 26) 

In the case (3 = 2 (unitary symmetry) the generating function ^• n, ((— s, s); £) can be ex- 
pressed Fredholm determinant 

Ef*i(-s,s);0 = det(l -£K ( ™; s) ) (3.27) 

where a \ is the integral operator on (— s,s) with kernel 

, , ( J a+l/2^x) J a _ 1/2 (Try) - J a+1/2 (7ry)J a _ 1/2 (7rx)) 
K™-(x,y) := (Trs) 1 / 3 ^)!/* A >- (3.28) 

evaluated at a = 1. Following the strategy which leds to ()3.22[) . the Fredholm determinant Q3.27JI 
for general a £ Z>o can be characterized as the solution of a nonlinear equation. Explicitly jllj 

E^i(-s, s);0 = exp ( dt ) (3.29) 

where a a satisfies the nonlinear equation 

{so'*) 2 + 4(-a 2 + sa' a - a a ){{a' a f -{a- (a 2 - sa' a + , a ) 1 / 2 } 2 ) = (3.30) 
subject to the boundary condition 

2(s/4) 2a+1 



s _o+ r(l/2 + a)r(3/2 + a)' 

In the case a = 0, (|3.28jl reduces to the sine kernel and the differential equation (|3.3Ujl reduces 
to l|3.4j) . For general a the differential equation (|3.3U|) is satisfied by an auxilary Hamiltonian 
for PHI (as distinct from PHI') [35]. 
Substituting (l3~2T)l) in (EOfl) gives 

a a (2iTs;0 [ 2 ™ a a (t;Q 
p 2 (s) = exp / at . (3.31) 



2ns 



a=£=l 



An application of this result can be made to the study of the zeros of the Riemann zeta function 
on the critical line (Riemann zeros). We recall that the Montgomery-Odlyzko law states that 
the statistics of the large Riemann zeros coincide with the statistics of bulk eigenvalues of an 
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nn (t) 



1.5 - 




t 



Figure 3: Comparison of nn(t) := p\ ' n '(s) for the matrix ensembles with unitary symmetry in 
the bulk (continuous curve) and for 10 6 consecutive Riemann zeros, starting near zero number 
1 (open circles), 10 6 (asterisks) and 10 20 (filled circles). 

ensemble of random matrices with unitary symmetry, where both the zeros and eigenvalues are 
assumed to be unfolded so as to have mean spacing unity. As a test of this law, in the 
empirical determination of P2' n '{s) for large sequences of Riemann zeros, starting at different 
positions along the critical line, was compared with (|3.31|) . The results, which are consistent 
with the Montgomery-Odlyzko law, are reproduced in Figure |21 A significant feature is that the 
empirical determination of P2' n (s) for the Riemann zeros is so accurate that it is not possible 
to compare against an approximate form of P2' n '(s) for the random matrices. Thus the exact, 
readily computable, Painleve evaluation (|3.31jl is of a practical importance. 
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4 Gap probabilities from the Okamoto r-function theory 

4.1 Other stragies 

The method of Tracy and Widom may be described as being based on function theoretic prop- 
erties of Fredholm determinants. Alternative methods which also lead to the characterization of 
gap probabilities in terms of the solution of nonlinear equations have been given by a number of 
authors. One alternative method is due to Adler and van Moerbeke [SHI; wno base their strat- 
egy on the fact that for suitable underlying weight W2, gap probabilities in the case of a unitary 
symmetry satisfy the KP hierarchy of partial differential equations known from soliton theory. 
The first member of this hierachy is then used in conjuction with a set of equations referred to as 
Virasora constraints, satisfied by the gap probabilities as a function of the endpoints of the gap 
free regions, to arrive at third order equations for some single interval gap probabilities. These 
third order equations are reduced to the u-form of the Painleve theory, making use of results 
of Cosgrove j31[2]. Borodin and Deift fl] have given a method based on the Riemann-Hilbert 
formulation of the resolvent kernel (|3.16|) |24j . This makes direct contact with the Schlesinger 
equations from the theory of the isomonodromic deformation of linear differential equations, and 
is thus closely related to the work of Jimbo et al. [23]. The other approach to be mentioned is due 
to Forrester and Witte [Tl|- It is based on Okamoto's development of the Hamiltonian approach 
to Painleve systems, and proceeds by inductively constructing sequences of multi-dimensional 
integral solutions of the a form of the Painleve equations, and identifying these solutions with 
gap probabilities for certain random matrix ensembles with unitary symmetry. 

For detailed accounts of all these methods, see 8 ( Ch. 6&7]. In the remainder of these 
lectures we will restrict ourselves to results from the work of Forrester and Witte which relate 
directly to gap probabilities in the bulk. 

4.2 Direct calculation of spacing distributions 

We have taken as our objective the exact evaluation of the bulk spacing distributions for the 
three symmetry classes of random matrices. So far exact evaluations have been presented not 
for the spacing distribution itself, but rather the corresponding gap probability, which is related 
to the spacing distribution by 1)2.2(1 . It was realized by Forrester and Witte ^S] that in all three 
cases one of the derivatives could be performed analytically by using theory relating to the a 
form of the Painleve transcendents. 

As an explicit example, consider the result (|ci.24|) . It was shown in ^5] that 




(4.1) 



where u satisfies the nonlinear equation 



s 2 {u") 2 = (4(u') 2 
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subject to the boundary condition 



8 s 2 8s 5 / 2 



u(s) ~ — h 



+ o+ 3 45 135vr 
Recalling now (|2.2j) we see that 



, ulk/n _ 2u((vr S /2) 2 ) _ ( u{t) 



? 

rf ulk (0; s) = Z^lljLZLJ. eX p (-J ^ dtj (4.2) 

(cf. (Z3Z» - 

The identity (|4.1j) can be understood from the approach to gap probabilities of Forrester and 
Witte. The key advance from earlier studies is that the generating function (|2.4|) . with p given 
by (|2,14|) . can be generalized to the quantity 

/co poo N 

dxi--- / dx N '[{( l -^X(l s /2^/2))\ s / 2 ~ a NXi\' J 'p(a N x 1 , a N x N ) 
-CO J — CO 

(4.3) 

and still be characterized as the solution of a nonlinear equation. This is also true at the hard 
and soft edges, and in the neighbourhood of a spectrum singularity (before the generalization 
the latter is controlled by the kernel (13.281) ). 

It is the generalization in the case of the hard edge which leads to (|4.1|) . The quantity of 
interest is defined by 

^«o, = Jim 7^y^((0, ^); f ; (« - s^>.) (") 

where 

N 

ly«i):= I cl.r i. • • • / f/.r.\ [ | e _Xi a;f f[ (** - Xj) 2 



r(o) := / dxi- - dx N Y[e~ Xl xf JJ (a* 

^° Z=l l<j<k<N 



L<j< 

(the factor I^f(a)/I^(a + /x), which is readily evaluated in terms of gamma functions, is chosen 
so that when s = 0, (|4.4|) is equal to unity). By using theory from the Okamoto r function 
approach to the Painleve systems PV and PHI' it is shown in |16j that 



E^ d ((0,s);^0 =exp / a, £) 

J o 

where u ft satisfies the differential equation 



dt 

7' 



(W) 2 - (a* + a) 2 (n') 2 - «'(4u' + l)(u - tv!) - ^±A U > _ !L = . (4.5) 

2 4^ 

Thus we have 

1 d / f s dt\ s a ( f s dt\ 

"e^ exp W a ' ^ 7) = 2 2 ^ 2 r(a + i)r(a + 2) exp ( J « A (*i«,M;OI^. T J, 

(4.6) 

which in the case o = — 1/2 reduces to (|4.1|l . 
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We also read off from Q4.6JI that 

- exp ^ u(t; a; j %(2 = -^s) 2 exp ( - ^ v(t)j) (4.7) 

where v(t) = — u h (t; a, /x; 1 0=1/2,^=2,^=1 an d thus satisfies (|4.5() appropriately specialized. The 
boundary condition consistent with (|4,7|) is 

gW ^« (1 + OW) + 55 ^L.(l + OW ). (4.8) 
Hence, according to Q2.34JI and l|2.2|) . 

^ ulk (0; S ) = 2pf lk (0;2 S ) + ^l^ns) 2 ) - l) exp ( - S(t)y). (4.9) 

The Okamoto r-function theory of PVI and PV allows (|4.3j) to be computed for general fi, 
and also its generalization in which there is a further factor | — s/2 — ajv2;| a in the product over 
I in the integrand 17 . These results allow not only the first derivative with respect to s of (|3.3j) 
to be computed by an identity analogous to (|4.1jl . but also the second derivative. In particular, 
it is found that 

P2 Ulk (0; s) = -js 2 exp J „(t; j (4.10) 

where v satisfies the nonlinear equation (which can be identified in terms of the a form of the 
PHI' equation) 

{sv"f + (v- sv'){v - sv' + 4 - 4(v') 2 } - W{v') 2 = 



subject to the boundary condition 



s^o 15 



The exact evaluations (|4.2|1 . (|4.9|) and (|4.1U|) are perhaps the most compact Painleve evalu- 
ations possible for the bulk spacing distributions. A striking feature of Q4.2JI and Q4.1UJI is that 
they are of the functional form a(s) exp(— J* Q S b(t) dt) and thus extend the Wigner surmise (|2.12j) 
and its (3 = 2 analogue in (|2.13|) to exact results. 
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